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ABSTRACT 

The Vela pulsar is one of the most exciting 7-ray sources and has been at the 
forefront of high-energy pulsar science since the detection of 7-ray pulsations at 
the radio period by SAS-2 in 1975. With the unprecedented angular resolution, 
effective area, field of view, and timing resolution, in the GeV band, of the Large 
Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope, the light curve 
of the Vela pulsar can be studied in greater detail than ever before. Using a 
timing solution derived solely from the LAT data, phase aligned with the radio 
emission, the spectrum of the Vela pulsar has been fit in intervals as small as 
0.0016 in phase. Significant variation is seen in the cutoff energy and photon index 
across the light curve, strongly supporting curvature radiation as the source of 
the high-energy 7-rays from the Vela pulsar. 

1. INTRODUCTION 
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The Vela pulsar has a long history 
of 7-ray observations dating back to the 
point source discovery and pulsed detec- 
tion with SAS-2 (Thompson et al.l974 
and Thompson et al.l975, respectively). 
Phase-resolved spectral re sults followed 



using data from COS-B (IGrenier et al. 



19881) and EGRET (Kanbach et al.l994 



and Fierro et al.l998) suggesting signifi- 
cant variation of the spectrum across the 
pulse, with the hardest emission observed 
between the two main 7-ray peaks. 

With the advent of the Fermi Gamma- 
ray Space Telescope {Fermi), it is now 
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possible to study the 7-ray emission from 
the Vela pulsar in unprecedented detail. 
The main instrument aboard Fermi is the 



Large Area Telescope (LAT) (lAtwood et al 



20091) ■ Using early calibration data (lAbdo et al 
2009al ) and ~2 months of survey data 
the LAT produced the highest resolution 
high-energy (HE, > 0.1 GeV) light curve 
of the Vela pulsar to date, identified a 
third peak which moves with increasing 
energy, and ruled out low-al titude mag- 



netosp heric emission models (lAbdo et al. 



2009bl ). With 11 months of survey data it 
is now possible to perform the finest phase- 
resolved study of the Vela pulsar to date by 
fitting the spectrum in 101 variable-width 
bins across the pulse. 

2. OBSERVATIONS 
2.1. Timing 

LAT data was used to construct a tim- 
ing solution for the Vela pulsar with 63 
us residua l s usiii g techniques described in 
Ray et al.l (120091 ). The 7-ray timing so- 
lutiorEI was phase aligned wit h data from 
the Parkes Radio Telescope (IManchester 
2008h . 



2.2. Data Selection 

LAT event data from 4 August 2008 to 
4 July 2009 with reconstructed energies 
from 0.02 to 100 G eV and belonging t o 
the "Diffuse" class flAtwood et al.i l2009h . 
as defined under the P6_V3 instrument re- 
sponse functions, were selected from within 
15° of the radio position. Time periods 



when the limb of the Earth infringed upon 
the region of interest and when the rock- 
ing angle of the LAT exceeded 43° were 
excluded to reduce background. Events 
were phase-folded using the Fermi plug-in 
now provided with the TEMPO^I software 
(IHobbs et al.l l2006l ). The resulting 7-ray 
light curve is shown in Figure 1, a detailed 
analysis of the ligh t curve can be found in 
Abdo et all (l2009ch . 



3. SPECTRAL ANALYSIS 
3.1. Phase- Averaged 

Events with energies > 0.1 GeV were 
used for spectral analysis. The Fermi 
science Tool^ (STs) v9rl5p2 were used 
to perform a binned maximum likelihood 
analysis (Cash 1979 and Mattox et al. 
1996) modeling all point sources found 
above the background with a test statistic 
> 25 in a preliminary versio n of the IFGL 
catalog (lAbdo et al.l l2009dl ): an extended 
source at the position of the Vela X pulsa r 
wind nebula (PWN) flAbdo et al.ll2009eh : 
and the v02 diffuse backgrounds, included 
with the Fermi STs. The phase-averaged 
spectrum. Figure 2, is best fit as an ex- 
ponentially cutoff power law. Equation 1, 
with the h parameter < 1: 



= No 

dE ° 



E \-r 



1 GeV 



exp 



E 



(1) 

likely due to the superposition of many 
spectral components with 6=1 and vary- 
ing values of E^; and F through the pulse. 
The best-fit spectral parameters are given 
in Table 1. 



*The timing solution will be made available 
through the Fermi Science Support Center 
[http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 



'l'http://tempo2. sourceforge.net/ 

•'■http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/ 
overview.html 
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Fig. 1. — Folded light curve of the Vela pulsar for events with reconstructed energies > 0.02 
GeV and within max[1.3°,1.6° — 3 * log]^Q(E/l GeV)] of the radio position, each bin has 750 
counts. 



3.2. Phase-Resolved 

The light curve was divided into 101 
variable-width bins with 1500 events each 
(> 0.1 GeV) using the energy- dependent 
angular selection described in the caption 
of Figure 1. The smallest bin has a width 
of 0.0016 in phase (~ 142 /is). The spec- 
trum of the Vela pulsar was fit assuming 
a simple exponentially cutoff power law, 
h = 1, in each phase bin. Normaliza- 
tions of the diffuse backgrounds and point 
sources within 5° of the pulsar were also 
kept free. The observed phase trends in 
r and E(7 are shown in Figures 3 and 4, 
respectively. 

Significant variation is seen in both 
parameters, confirming that the hardest 
emission is between the two peaks. Ec 
rises sharply through the main peaks and, 
suprisingly, between the peaks as well, fol- 
lowing the position of the third peak ob- 



served at higher energies. Phase-resolved 
analyses of EGRET data suggested a dras- 
tic change in F through the two main 
peaks; however, the LAT data does not 
confirm this, instead finding F to be very 
consistent with a constant value through 
both peaks. 

To evaluate the significance of the fea- 
tures in Figures 3 and 4, the pulsar and 
surrounding region were simulated using 
the Fermi ST qto hssim and the built- 
in Pulsar Spectrum (IRazzano et al.l l2009l ). 
The LAT-only timing parameters and 
spectral parameters from a phase-averaged 
fit with b = 1 were used as input to the 
simulation. The Vela X PWN was included 
as a point source. The simulation suggests 
that point-to-point variations of 0.6 GeV 
in Ec and 0.05 in F should be expected 
from the fitting technique. As such, point- 
to-point variations less than these values 
can not be considered significant. 
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Table 1 

Vela Phase- Averaged Spectral Parameters'' 



No (10-6 cm"2 GeV"') 


3.73 ±0.31 ±1.04 


r 


1.37 ±0.03+^:^^ 


Ec (GeV) 


1.31 ±0.18+^:^ 


h 


0.68±0.03+^:i^ 


Flux (0.1-100 GeV) (10"^ cm-^ g-^) 


1.07 ±0.01 ±0.03 



'^First errors are statistical, second are systematic. 




1 10 
Energy (GeV) 

Fig. 2. — Phase- averaged 7-ray spectrum 
of the Vela pulsar. Best-fit model with h < 
1. The representative data points are from 
likelihood fits to each energy range with 
the pulsar spectrum modeled as a power 
law. 
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Fig. 3. — Photon index vs. phase. Errors 
are statistical only. Results are only shown 
for bins in which the pulsar was found abve 
the background with a test statistic > 25. 

maxima near the start of the trailing edges. 



To better evaluate the behavior of Ec 
in the two main peaks the fit was repeated 
but with P fixed to the best-fit values of 
1.72 ±0.01 and 1.58 ±0.01 for the first and 
second peaks, respectively. The resulting 
trends in Ec are shown in Figures 5 and 6 
for the first and second peaks, respectively. 
Ec rises fairly smoothly in both peaks with 



4. Discussion 

The HE 7-ray emission from pulsars 
has been theorized to be curvature radi- 
ation (CR) from electrons (or positrons) 
accelerated along the magnetic field lines 
by the parallel component of the elec- 
tric field (E||) (Romani 1996 and Hi- 
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Fig. 4. — Cutoff energy vs. phase. Er- 
rors are statistical only. Results are only 
shown for bins in which the pulsar was 
found above the background with a test 
statistic > 25. 



Fig. 5. — Cutoff energy vs phase for the 
first peak. The photon index was held fixed 
to the best fit value of 1.72 ± 0.01 between 
the dashed vertical lines which correspond 
to the phase range 0.112 < (p < 0.155. 



rotani & Shibata 1999). The electrons 
will be CR reaction limited such that a 
steady-state Lorentz factor is maintained. 
In outer-magnetospheric emission models 
(e.g. Muslimov & Harding 2004, slot gap 
(SG); Zhang et al. 2004; and Hirotani 
2008, outer gap (OG)) Ey depends on the 
value of the magnetic field at the light 
cylinder and the gap width. These mod- 
els give similar cutoff energies. Equation 2 
(units of mc^), ranging from 1 to 5 GeV, 
consistent with what is observed in Vela 
and other 7-ray pulsars. 

EcR = --7Sr = 0.32Ae( J) 'pi (2) 
2 pc V e / 

Interestingly, the cutoff energy (Ecr) 
depends on the local field line radius of 
curvature (pc)- Emission across the pulse 
originates from different ranges of emission 



radii (Romani 1996 and Cheng et al. 2000 
(OG) and two-pole caustic (TPC) Dyks & 
Rudak 2003), implying that phase- resolved 
spectroscopy should map out the emission 
altitude. Large variations of pc with phase 
are also expected in these models and map- 
ping the minimum pc, using basic geomet- 
ric models, can produce trends similar to 
what is seen in Figure 4, but full radia- 
tion models will be needed to match all of 
the features. As the LAT continues to ac- 
cumulate more events from Vela it will be 
possible to map the HE 7-ray emission re- 
gions in more detail and better understand 
the radiative processes involved. 
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